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Abstract 

Existing  models  to  obtain  intensity  and  phase  changes  of  sound  due  to  thermal 
fluctuations  in  the  ocean,  such  as  that  developed  by  Chernov,  are  limited  in  their 
application.  They  are  valid  only  for  an  isotropic  medium  and  where  the  range  is  much 
greater  than  the  correlation  distance  of  the  fluctuations.  Extensions  to  short  range  and  for 
an  anisotropic  medium  are  described  here,  both  for  the  stochastic  wave  and  ray  models. 
Theoretical  expressions  obtained  for  acoustic  fluctuations  are  easily  evaluated.  Numerical 
estimates  are  provided  using  stochastic  wave  arid  ray  models  for  certain  parameters  of 
range,  frequency  and  temperature  variation,  of  interest  in  target  strength  measurements. 
These  estimates  indicate  acoustic  fluctuations  for  an  anisotropic  medium  are  always  higher 
than  for  an  isotropic  media.  For  a  range  of  100  m,  frequency  of  1 00  kHz  and  temperature 
fluctuation  of?  xlO~^  C  deg,  the  root  mean  square  relative  intensities  and  the  root  mean 
square  phase  fluctuations  are  on  the  average  higher  by  a  factor  of  6,  where  the  vertical 
correlation  lengths  are  the  same  and  the  vadial  correlation  length  is  about  50  times  larger. 
Closed  form  expressions  in  the  deterministic  ray  approximation  are  given  to  evaluate 
intensity  and  phase  changes  due  to  the  mean  temperature  gradient.  The  validity  of  the 
deterministic  ray  theory  results  is  discussed.  Results  from  the  stochastic  wave  and  ray 
models  are  examined  for  the  applicability  of  the  stochastic  ray  model  at  high  frequencies. 

At  mega  hertz  frequencies  and  ranges  of  1-3  metres,  which  are  of  interest  in  acoustic  mine 
imaging,  estimates  of  the  phase  fluctuations  from  the  stochastic  ray  approximation  and  the 
wave  model  are  indistinguishable.  Theoretical  expressions  of  the  phase  structure  function 
are  giiren  in  the  ray  approximation  to  provide  estimates  of  the  phase  coherence  between 
sigtuds.  These  results  can  provide  valuable  comparisons  with  experimental  results,  when 
such  measuremetUs  are  carried  out  unth  frequencies  in  the  mega  hertz  band  and  at  ranges 
of  1-3  metres. 
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Models  to  Estimate  High  Frequency 
Acoustic  Scattering  Due  to  Thermal  Fine- 
and  Micro-Structure  of  the  Ocean 

1.  Introduction 

Scattering  in  the  ocean  can  result  from  several  sources.  These  include  organic  and 
inorganic  particles  (such  as  plankton  and  sand),  air  bubbles,  the  sea  bottom  or 
surface,  and  chang*  s  in  current  or  temperature.  This  paper  focuses  attention  on 
high  frequency  acoustic  scattering  due  to  temperature  changes  in  the  ocean. 

Scattering  of  sound  due  to  temperature  changes  can  cause  attenuation  and 
fluctuations  in  amplitude  and  phase.  Amplitude  fluctuations  cause  erratic 
detection  ot  targets,  since  the  targets  tend  to  be  detected  when  the  signal  is  strong 
and  lost  when  the  signal  is  weak  [Urick,  1983,  pl83].  Further,  in  acoustic 
underwater  imaging,  temperature  fluctuations  can  cause  varying  intensities  and 
phase  errors  at  different  positions  of  the  image,  and  loss  of  correlation  between 
signals.  Thus  the  iiUportance  of  this  study. 

In  this  paper  a  number  of  questions  are  looked  at,  in  relation  to  high  frequency 
acc»ustic  scattering  due  to  temperature  changes  in  the  ocean.  How  do  we  obtain 
reliable  theoretical  estimates  of  acoustic  attenuation,  intensity  and  phase 
fluctuations  from  ocean  temperature  measurements?  Can  we  theoretically 
estimate  intensity  and  phase  fluctuations  due  to  temperature  changes  in  target 
strength  measurements  where  the  range  is  around  100  m?  In  the  study  of 
underwater  acoustic  imaging  the  range  involved  is  a  few  metres.  Is  it  possible  to 
obtain  theoretical  estimates  at  this  short  range?  The  ocean  being  an  anisotropic 
medium,  how  does  one  obtain  estimates  of  acoustic  fluctuations,  since  most  of  the 
models  treat  it  as  isotropic.  How  good  is  the  stochastic  ray  approximation  at 
frequencies  of  100  kHz  and  above?  An  experiment  usually  carried  out  is  the 
measurement  of  the  coherence  of  the  phases  of  an  acoustic  signal  at  two  receivers 
separated  in  space.  The  phase  structure  function  is  a  measure  of  this  coherence.  Is 
it  possible  to  obtain  theoretical  estimates  of  this  quantity  at  short  range? 

Flatte  [1979,  p931  presented  the  different  sound  transmission  regions  in  the 
A-^space.  Here  A  is  the  diffraction  parameter,  and  <I>  the  root  mean  square 
(rms)  phase  fluctuations  in  the  ray  model.  These  parameters  are  defined  later  in 
this  report.  Fig.l  is  a  diagram  of  this  space.  The  unsaturated  region  corresponds  to 
weak  scattering,  while  the  saturated  region  is  that  of  strong  scattering.  The  ray 


region  is  where  diffraction  effects  are  small.  For  an  isotropic  medium  expressions 
for  A  [Flatte,  1979,  p91,  eqn.6.2.3]  and  <I>"  [Flatte,  1979,  p92,  eqn  6.2.5]  are 
available.  The  expression  for  <!>'  is  valid  only  when  the  range  (L)  is  very  much 
greater  than  the  correlation  length  ). 


Figure  1;  A  -  <I>  space,  where  A  is  the  diffraction  parameter,  and  ih  is  the  rms  phase 
fluctuation  in  the  ray  model. 


In  the  ray  region,  expressions  for  the  relative  intensity  fluctuations  and  phase 
fluctuatiorts  are  available  orUy  for  an  isotropic  medium,  in  a  domain  where  L  is 
very  much  greater  than  [Chernov,  1960,  p34,  p29;  Flatte,1979,  p92].  No  closed 
form  results  tue  available  for  short  range  where  L  is  less  than  or  equal  to  and 
for  an  anisotropic  medium.  Again  in  the  unsaturated  Rytov  region,  expressions 
for  relative  intensity  and  phase  fluctuations  are  available  only  for  L  very  much 
greater  than  and  for  large  [Chernov,  1960,  p83].  No  closed  form 
expressions  are  available  for  L  less  than  or  equal  to  ,  or  for  an  anisotropic 
medium. 

The  paper  discusses  extensions  to  these  models  to  obtain  estimates  of 
fluctuations  at  all  ranges  in  the  ray  and  Rytov  regions,  both  for  isotropic  and 
anisotropic  media.  Numerical  estimates  are  provided  for  some  situations  to  give  a 
measure  of  the  magnitude  of  these  fluctuations.  Comparisons  are  made  between 
the  Chernov  results  md  the  results  of  the  stochastic  ray  model.  Mathematical 
detail  is  kept  to  a  minimum;  the  emphasis  will  be  on  the  application  of  models  to 
obtain  acoustic  estimates. 

The  structure  of  the  paper  is  as  follows.  Some  experimental  temperature  and 
temperature  gradient  profiles  are  examined.  Using  deterministic  ray  theory, 
closed  form  results  for  intensity  and  phase  changes,  with  reference  to  a  medium 
where  there  is  no  temperature  gradient,  are  given.  These  results  provide  a 
qualitative  imderstanding  of  the  effect  of  the  mean  temperature  gradient  on 
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acoustic  transmission,  and  quantitative  estimates  of  intensity  and  phase  changes 
under  certain  limited  conditions.  Next,  a  statistical  model  based  on  wave  theory  is 
developed  and  examined.  Theoretical  expressions  for  the  mean  square  scattered 
pressure,  attenuation,  intensity  and  phase  variance  at  a  "•t  separated 
horizontally  from  the  transmitter  are  given.  Comparisons  are  made  with  other 
models  suclt  as  the  Chernov  wave  model  and  a  new  stochastic  ray  model.  Finally 
a  theoreticad  expression  for  the  phase  structure  function  based  on  the  stochastic 
ray  model  is  described. 

Velocity  fluctuations  that  occur  in  the  ocean  can  be  classified  in  general  into 
three  time  scales  (Anderson  and  Zahuranec,  1977,  pl35J.  Velocity  fluctuations  that 
arise  due  to  earth's  rotation  occur  in  time  scales  in  excess  of  1  /  /  hours,  where  / 
is  the  inertial  frequency  or  Coriolis  frequency  in  /trs*' .  1  /  is  approximately  24 
hours.  Velocity  fluctuations  that  arise  due  to  density  gradient,  stratification,  or 
internal  waves  occur  in  time  scales  between  1  /  to\ !  ,  where  corresponds 

to  the  stability  frequency  or  the  Brunt  Vaisala  frequency  ( 1 !  j\,  lies  approximately 
between  30  mins  to  24  hrs).  Velocity  fluctuations  due  to  temperature  fine  structure 
and  micro  structure  occur  typically  in  a  time  scale  less  than  30  mins.  The  length 
scale  of  the  structure  .is  the  basis  of  division  between  fine  structure  and  micro 
structure.  Fine  structure  corresponds  to  a  length  scale  between  1  m  and  100  m 
while  for  the  micro  structure  it  is  below  1  m  [Flatte,  1979,  pl5]. 


2.  Some  Experimental  Temperature  and 
Temperature  Gradient  Profiles 

Fig.  2(a)  is  an  example  of  a  temperature  profile  obtained  off  Perth,  Western 
Australia  (from  the  RV  Franklin).  Fig.  2(b)  is  the  corresponding  temperature 
gradient  profile.  Temperature  is  given  at  every  2  m  interval  and  has  an  accuracy 
of  10*’  C  deg.  The  plots  indicate  that  the  temperature  gradient  has  more  structure 
than  the  temperature  prc  ile.  In  the  next  section  it  will  be  shown  that  it  is  the 
temperature  gradient  that  is  responsible  for  acoustic  refraction. 

Fig.  3(a)  is  an  example  of  the  temperature  profile  obtained  at  Woronora  dam,  a 
fresh  water  lake  approximately  50  km  south  of  Sydney,  New  South  Wales.  Fig. 
3(b)  is  the  temperature  gradient  of  this  profile.  The  profile  is  over  a  depth  interval 
of  10  m.  One  can  observe  here  some  temperature  inversions  that  are  characteristic 
of  surface  layers.  Normally,  temperature  decreases  with  depth  {dT /  dZ  <0). 

Fig.  4  shows  the  temperature  variation  with  time  over  the  same  10  m  depth 
interval  at  three  times.  Plotted  here  is  the  temperature  difference,  that  is  the 
temperature  at  time  t  minus  the  temperature  at  t=v),  where  the  temperature  at  t=0 
corresponds  to  the  profile  shown  in  Fig.  3(a).  The  three  parts  of  Fig.  4  correspond 
to  variations  after  3, 6  and  9  mins  respectively.  The  diagram  shows  that  even  at 
the  small  time  scales  considered  here  there  is  a  temperature  variation,  the 
magnitude  of  which  on  the  average  varies  from  10‘‘  to  10''c  deg  per  minute. 
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Figure  3:  (a)  Temperature,  (b)  Temperature  gradient  profile  obtained  at  Woronora  dam, 
approximately  50  km  south  of  Sydney,  New  South  Wales. 
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Figure  5;  Discrete  probability  fuiictwn  of  the  temperature  gradient  obtained  off  Perth  in 
Western  Australia,  and  A  times  the  gamma  probability  density  function  ( the  smooth 
curve).  A  =  10'-^. 


Fig.  5  is  the  discrete  probability  function  of  the  temperature  gradient  obtained 
using  the  Perth  data  shown  in  Fig. 2.  This  was  obtained  from  the  frequency  of 
occurrence  of  a  particular  temperature  gradient  (using  the  Numerical  Algorithm 
Group  (NAG)  routine  GOl  AEFy  For  ease  of  comparison  with  the  model  of  Hayes 
et  al  (1975),  the  temperature  gradient  plotted  is  the  negative  of  the  observed 
temperature  gradient.  The  probability  .'unction  shows  a  large  skewness  and  can  be 
fitted  quite  well  to  a  function  of  the  form  A  /(x),  where  A  is  the  discretisation 

— JC 

interval  (equal  to  10*’  in  this  example)  and  f(x)  =  which  is  the 

gamma  probability  density  function.  The  variable  X  corresponds  to  the 
temperature  gradient.  The  mean  and  vanance  of  the  X  distribution  give  the 
constants  and  Pq  [Freund,  1965,  pl47I,  while  the  result  that  the  integral  of  the 
ganuna  probability  density  function  over  all  values  of  X  equals  unity  gives  the 
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constant  .  The  model  proposed  by  Hayes  et  al  [1975]  is  the  basis  of  the  choice 
of  f{x).  This  model  assumes  layers  of  water  that  are  nearly  homogeneous, 
separated  by  sheets  of  sharp  gradient  with  the  mixing  process  occurring 
randomly  in  the  vertical  according  to  Poisson’s  law.  This  model  is  the  basis  of  the 
choice  of  the  anisotropic  correlation  function  used  later  in  this  paper.  So  far,  only 
vertical  temperature  profiles  have  been  presented.  However  to  obtain  accurate 
estimates  one  also  requires  temperature  profiles  at  horizontal  separations. 

Next,  the  effect  of  the  mean  temperature  gradient  on  the  acoustic  intensity  and 
phase  changes  will  be  examined,  using  simple  deterministic  ray  theory.  The 
reference  is  a  medium  where  there  is  no  temperature  gradient. 


3.  Deterministic  Ray  Theory  Results 

The  speed  of  sound  (C  )in  the  ocean  can  be  given  by  an  empirical  equation  of  the 
form, 

c  =1449.2  +  4.67-0. 0557^  +0.000.37^  +(1.34  -  0.0127)(S-.35)  +  0.016Z 


where  T  is  the  temperature  in  degrees  Centigrade,  S  is  the  salinity  in  parts  per 
thousand  and  Z  the  depth  in  metres  [Medwin,  1975].  The  speed  of  sound  in  the 
ocean  is  therefore  a  function  of  temperature,  salinity  and  depth.  However,  as  seen 
from  the  above  equation  the  largest  coefficient  is  that  of  temperature,  which 
indicates  that  temperature  change  is  the  dominant  cause  of  velocity  variation. 

Consider  the  case  where  we  can  approximate  the  above  equation  by  c  =  Cq  +  ^2 


where  g  is  equal  to  the  velocity  gradient 


dc 

7z 


.  This  g  is  mainly  due  to 


temperature  variation.  From  ray  theory,  the  path  of  the  ray  is  given  by  [Officer, 
1958,  p48]; 


ds  ^  dZ 


where  p  = - ^,s  is  the  ray  path,  0y  ti\e  angle  at  the  source,  and  the  speed  of 

sound  near  the  source.  We  measure  6  and  as  the  angles  of  the  ray  made  with 

d.'T 

the  vertical  as  shown  in  Fig.  6(a),  with  depth  increasing  downwards.  If  — 

changes  with  depth,  then  —  changes  with  depth. 

dZ 


If  — >0 

dZ 


then 


‘<9  n 

ds 


that  is  0  increases  with  S  .  In  other  words  the  ray  bends  upwards.  On  the  other 
hand. 


if 


dc 

dZ 


<0 


then  —  <0, 
ds 


that  is  6  decreases  with  5 ,  and  the  ray  bends  downwards.  Thus  due  to 
temperature  gradient,  and  therefore  velocity  gradient,  rays  undergo  refraction. 
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Let  us  now  look  at  intensity  changes  due  to  refraction.  Suppose  we  have  a 
receiver  (/^ )  and  transmitter  ( T)  at  the  same  depth  but  separated  horizontally  by 
a  distance  L,  as  indicated  in  Fig.  6(a).  Let  us  assume  always  that  the  transmitter  is 
at  z  =  0  and  the  speed  of  sound  near  it  is  .  For  a  medium  with  a  linear  velocity 
gradient  g,  the  intensity  loss  due  to  refraction  measured  in  decibels  (dB)  is  given 
by 

A/(Ri)(dB)  =  10Lo^lo(l  +  ^^).  (1) 

4co^ 


This  was  obtained  from  Officer's  equations  (2-83),  (2-85)  and  from  the  expression 
for  the  intensity  at  horizontal  range  L  given  by. 


1  = 


Eq  sin  00 
Lcose(^) 


(2) 


[Officer,  1958,  p50,  eqn.2-54].  In  eqn(2),  E^^  is  the  amount  of  energy  emitted  per 
unit  time  per  unit  solid  angle.  These  equations  are  invalid  for  either  L  equal  to  0 
or  for  0Q  equal  to  0,  as  seen  from  Officer's  eqns  (2-54),  (2-83)  and  (2-84).  Equation 
(2)  also  fails  at  0  =  0  5n,  that  is  at  a  turning  point.  The  intensity  loss  A/(/?, )  is 
the  difference  in  intensity  at  ,  between  a  medium  with  a  constant  velocity 
gradient  g  and  one  where  g=0.  The  larger  the  magnitude  of  g,  the  bigger  is  the 
intensity  change.  Another  feature  to  be  noted  is  that  the  intensity  changes  due  to 
refraction,  when  receiver  and  transmitter  are  at  the  same  horizontal  distance,  is 
independent  of  the  sign  of  g'.  What  this  means  is  that  although  the  direction  of 
bending  is  different,  this  being  dependent  on  the  sign  of  'g',  the  amount  of 
bending  at  ,  that  is  the  normal  cross  sectional  area  at  7?, ,  is  the  same.  To  get  an 
estimate  for  the  intensity  losses  due  to  refraction,  let  us  look  at  some  results  for 
L  =  100  m  and  Cq  =  1500  m/s.  The  intensity  loss  for  g  equal  to  -0.23, 

-0.46  ^''are  2.55X10"^  dBand  1.02X10”^dB  respectively.  These  results  are 
obtained  using  eqn  (1).  The  velocity  gradients  correspond  to  temperature 
gradients  of  -0.05  C  deg/m  and  -0.1  C  deg/m  respectively.  Such  gradients  are 
observable,  as  may  be  seen  from  Fig.  2(b). 

Consider  the  case  where  the  transmitter  and  receiver  are  separated  by  a  distance 
R  and  at  the  same  time  they  are  displaced  vertically  by  Z .  We  again  take  a 
coordinate  system  where  the  depth  increases  downwards,  so  that  if  the  receiver  is 
above  the  transmitter,  z  <0,  while  if  the  receiver  is  below  the  transmitter,  z  >  0. 
The  intensity  loss  due  to  refractioti  in  this  case  is  given  by, 

z  z 

A/(R2)(rfB)  =  10log,o(l  +  ^-^  +  — )  (3) 

4c^  Co 


When  z  =  0  and  R  =  L,  eqn  (3)  reduces  to  (1).  The  result  given  in  eqn  (3)  is 
obtained  using  Officer’s,  eqn  (2-89),  (2-83)  and  (2-85).  Both  eqns  (1)  and  (3)  are 
invalid  either  for  L  =  0  or  for  0g  =  0 .  Unlike  eqn  (1),  for  z  0  the  intensity  loss 
depends  on  the  sign  of  g.  For  ^  <  0  (velocity  decreases  with  depth),  intensity  loss 
due  to  refraction  will  decrease  with  z  for  constant  R.  On  the  other  hand  for 
g>0  (velocity  increases  with  depth),  intensity  loss  due  to  refraction  will  increase 
with  Z  for  constant  R. 
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An  interesting  corollary  of  eqn  (3)  is  that  one  can  find  a  depth  Z  for  a  particular 
g  and  R  where  the  interrsity  change  due  to  refraction  becomes  zero.  In  other 
words  the  normal  area  of  the  ray  tube  is  the  same  as  that  for  an  isospeed  medium. 
This  depth  is  given  by 


z 


4^0 


For  g  >  0,  this  is  above  the  transmitter,  whereas  for  g  <  0  it  is  below  the 
transirutter. 

Let  us  now  consider  the  case  of  normal  incidence  where  6^=0.  We  begin  with 
Officers  eqn(2-81). 


ds 


sin  00 
^0 


g 


(4) 


For  0Q  =  0 ,  the  above  expression  becomes  zero,  which  implies  there  is  no 
refraction,  and  hence  no  intensity  or  phase  changes  due  to  refraction. 

—2c 

When  0Q  0,  but  Z.  =  0 ,  Officer’s  eqn  (2-83)  and  (2-85)  gives  z  = - 2. .  in  other 

s 

words  there  is  a  unique  position  where  the  ray  crosses  the  Z  axis.  At  this  position, 
using  Officer's  eqn  (2-80),  the  acoustic  velocity  c  =  — c^, .  Thus  the  situation  with 
00  ^  0 ,  but  L  =  0  is  not  realistic. 

Let  us  now  look  at  the  effect  of  the  mean  temperature  gradient  on  the  phase 
change,  with  reference  to  a  medium  where  there  is  no  temperature  gradient. 
Consider  the  case  where  the  transmitter  and  the  receiver  are  at  the  same  depth 
but  separated  horizontally  by  Z, .  In  this  case  derivation  of  a  closed  form  result  is 
possible.  Starting  from  Officer's  eqn  (2-84),  and  using  eqn  (2-83),  with  the 
trigonometric  relationship  between  the  full  and  half  angle  for  cotangent,  one  can 
show  that  the  phase  change  due  to  refraction  is 


A<>(Ri)  =  (27t/)  jln 


2c, 


4c" 


g*0 


where  /  is  frequency  in  Hz. 

Consider  g  =  -g^ .  Substituting  for  g,  one  obtains  after  some  algebra: 


I 


2 

In 

ft*' -'>^1 

2co  4Cp 

Co 

gi>0  (5) 


This  shows  that  the  phase  change  is  independent  of  the  sign  of  g.  For  small  (gj, 

L^lgl 

A«(R,)»(27tZ)— 

4co‘‘ 

showing  that  the  larger  the  value  of  |g|,  the  bigger  is  the  phase  change.  If  |g|  is 
equal  to  0,  then  there  is  of  course  no  phase  change  due  to  refraction. 
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Consider  the  case  where  the  transmitter  and  receiver  are  at  a  separation  R  but 
displaced  vertically  by  Z .  The  method  to  obtain  phase  change  due  to  refraction  is 
as  fic^ws.  For  a  particular  R  and  2  one  calculates. 


cx)t0o  = 


R^g  +  lzcp 
2co(R^-z^)^ 


This  was  obtained  using  Officer's  eqns  (2-51)  and  (2-83).  Starting  horn  Officer's 
eqn  (2-84),  and  again  expressing  the  cotairgeitt  of  the  half  angles  in  terms  of  the 
full  arrgles,  the  phase  change  is  given  by. 


6j^(R2)  =  (2nf) 


■ 

iln 

cotOo  +(cot^  00  +1)^ 

_R 

$ 

Ycot0+(cot^  0+1)^ 

Co 

(6) 


where  from  Officer's  eqn  (2-51), 


cot0  =  — 


c;  Co 


*^0 


Equations  (5)  and(6)  are  valid  for  g  0.  If  g  =  0  there  is  no  refraction  and  hence 
no  phase  change.  The  parameter  y  is  defined  as  follows: 


If 


{R^  -  >^cot{Qo) 
i 


then  7  =  -l  if  g>0  while  y  =  l  if  g<0. 

If  (R2-2^)^2<^cot(0o) 

$ 

then  Y  =  1  if  g>0  while  y  =  -l  if  ^<0. 
If 

S 

then  7  =  1,  and  ® 


For  2  =  0,  orie  can  show  that  eqn  (6)  reduces  to  eqn  (5).  For  2  ^  0  the  phase 
change  is  dependent  on  the  sign  of  To  uitderstand  this,  let  us  look  at  the  ray 
pathsdrawninFig.  6(b).  When  the  receiver  is  in  a  regitm  where  the  velocity 

is  greater  than  .  This  would  therefore  shorten  the  time  taken  by  die  ray  to  reach 
the  reemver,  whm  compared  to  a  medium  where  g  =  0.  This  results  in  a  negative 
phase  change.  As  depth  iitcreases,  so  does  the  phase  change,  where  the  reference 
is  a  medium  having  a  ocoistant  sound  speed  Cq.  If  g<0,  the  velocity  of  sound 
near  the  receiver  is  less  than  that  near  the  transmitter.  This  would  increase  the 
time  taken  to  reach  the  receiver,  when  compared  to  a  medium  where  g  =  0.  This 


IS 


would  result  in  a  positive  phase  change.  Increasing  the  depth  increases  the  phase 
change. 

The  only  discontinuity  and  singularity  that  occur  in  eqn  (6)  are  at  8q  =  0  and 
0  =  0  respectively.  If  0q  =  0 ,  it  has  been  shown  before  that  there  is  no  refraction, 
and  h«ice  no  phase  change.  If  0  =  0,  but  0Q  ^  0,  then  from  Officer's  eqn  (2-85), 

z  =  — ^  This  position  corresponds  to  c  =  0,  which  is  not  realistic.  In  other  words, 

if  0g  ^  0 ,  then  0  can  never  equal  zero. 

To  get  an  estimate  of  the  phase  changes  and  intensity  changes  due  to  refraction 
let  us  look  at  the  system  depicted  in  Fig.  6(c).  The  system  consists  of  three 
receivers,  all  at  the  same  horizontal  distance  L,  but  vertically  displaced  by  — z,0, 
and  Z  m.  In  imaging  it  is  not  intensity  or  phase  change  at  one  point  that  is 
important  but  the  differences  between  them.  The  phase  and  intensity  changes  at 
each  of  the  receivers  and  R,  are  calculated  using  eqns(l),  (3),  (5)  and  (6). 

Let  us  consider  three  cases.  The  first  case  corresponds  to  Z-  =  3  m,  2  =  0. 5  m, 
f  =  1  MHz  and  g  ~  -0. 23 .  The  intensity  changes  at  and  R^  are 
3. 33X10"^,  2. 29X10“^,  and  -3. 33X10~^  dB  respectively,  while  the  phase  changes 
are  -0. 488,  -1. 1 1X10“*  and  0. 488  radians  respectively.  The  next  case  considered 
is  Z,  =  10  m,  2  =  0. 5  m,  /  =  100  kHz  and  g  =  —0. 23 .  The  intensity  changes  at 

and  are  3. 36X10“*, 2. 55X10-^  and  -3.30A'10“‘  dB  respectively,  while 
the  phase  changes  are  -0.161,-4.10X10“*,  and  0.161  radians  respectively.  For 
the  third  case,  increasing  Z,  to  100  m  at  100  kHz  gives  intensity  changes  at 
R^,R^  and  R^  of  5.88X10^,2.55X10“*  and  -7.77X10“®  dB  respectively,  while 
the  phase  changes  are  -2.02,-0.41  and  1.195  radians  respectively. 

The  results  described  so  far,  for  intensity  and  phase  changes  using  ray  theory, 
were  obtained  by  considering  the  medium  to  be  one  where  the  speed  of  sound 
varies  linearly  with  a  constant  gradient  and  is  homogeneous  horizontally.  In  a 
medium  that  consists  of  many  vertical  layers  with  different  constant  velocity 
gradients,  it  is  impossible  to  obtain  closed  form  expressions  for  intensity  and 
phase  changes.  In  the  case  where  one  can  identify  the  presence  of  more  than  one 
layer,  the  solution  is  obtained  by  solving  numerically  Officer’s  eqns  (2-83  to  2-85) 
for  each  layer.  The  ray  solution  is  an  exact  solution  of  the  wave  equation  at  infinite 
frequency.  At  finite  frequencies  the  ray  solution  is  a  good  approximation  to  the 
wave  equation  if  (Tolstoy  and  Clay,  1966,  p52  and  Officer,  1958,  p391 


ko^A 


«1 


where  *b(=^%)  is  the  wave  number,  A  the  amplitude  and  the  Laplacian. 

The  ray  solution  fails  at  caustics  and  shadow  zones  [Tolstoy  and  Clay,  1966,  p57]. 
In  a  plane  stratified  medium,  like  the  one  assumed  here,  it  also  fails  at  turning 
points.  Eqits  (1),  (3),  (5),  (6)  are  obtained  under  the  assumption  that  there  are  no 
reflections  or  scattering  from  the  boundaries,  and  no  interference.  Under  such 
circumstances  for  a  plane  stratified  medium  the  equation  of  a  caustic  is  defined  by 

=  0  (Kravtsov  and  Orlov,  1990,  pl37],  where  L  is  the  horizontal 

range.That  would  be  the  equation  of  the  caustic  is  also  seen  from  eqn  (2).  Using 
Officer's  eqns  (2-83, 2-85)  for  L  aitd  differentiating  with  respect  to  0o,  the  equation 
for  the  caustic  can  be  obtained.  The  equation  obtained  is  an  implicit  equation  in 
Og.  For  a  single  layer  having  a  constant  velocity  gradient  (as  assumed  here). 
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numerical  solution  indicates  absoice  of  caustics.  This  can  also  be  shown 
analytically  as  follows.  For  a  medium  with  constant  velocity  gradient  the  ray 
paths  for  different  for  different  6^  are  circular  arcs  that  pass  through  the  source 

2ca 

widi  a  horizontal  range  of  — ^  cot  do  •  Since  the  horizontal  range  decreases 

monotonkally  with  Gg,  ray  paths  with  different  6^  do  not  intersect. 

For  near  horizontal  rays,  using  Officer’s  eqns  (2-27,2-46,2-83,2-85,2-13)  and  the 

assumption  that  and  ^  are  approximately  zero,  one  can  show  that  the  ray 
A  dz  aL 

solution  will  be  a  good  approximation  to  the  wave  equation  if 


„  ^  ,  .  1  3A  38  .  30  L  II 

For  steep  rays  ( where  z  is  large),  assuming  and  to  be  ail 

A  oL  oL  dz 

approximately  zero,  one  can  show  using  Officer's  eqns  (2-27, 2-46)  that  the  ray 
solution  will  be  a  good  approximation  to  the  wave  equation  U/»g  (Tolstoy 
and  Clay,  1966,  p50]. 

So  far  in  this  section,  expressions  to  obtain  intensity  and  phase  changes  in  a 
medium  with  a  mean  temperature  gradient  have  been  described.  The  reference 
has  been  a  medium  where  there  is  no  temperature  gradient.  The  usefulness  of 
these  simple  dosed  form  expressions  is  that  one  can  estimate  changes  in  intensity 
and  phase,  due  to  changes  in  the  mean  temperature  gradient,  in  regions  where  the 
ray  solution  is  valid. 


4.  Stochastic  Wave  and  Ray  Theory  in  a 
Random  Medium 

In  most  cases  in  the  real  ocean  one  can  write 

c  =  CQ+gz  +  small  random  part 

The  following  sections  look  at  obtaining  acoustic  estimates  such  as  intensity  and 
phase  fluctuations,  and  attenuation  from  the  random  part  of  the  sound  velocity. 
This  is  done  using  both  stochastic  wave  and  ray  models.  The  stochastic  ray  model 
suffers  from  many  deGciendes,  in  that  there  is  no  frequency  dependence  for 
intensity  loss,  and  ignores  diffraction  and  interference.  However  it  is  an  exact 
solution  at  infinite  frequency.  By  comparing  the  results  of  the  stochastic  ray 
model,  with  that  of  the  stochastic  wave  model,  the  region  of  applicability  of  the 
stochastic  ray  model  for  finite  frequencies  is  examined. 
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4.1  Correlation  Function,  Mean  Square  Scattered  Pressure  and 
Attenuation 


The  basis  of  this  model  is  the  approach  developed  by  Liu  [1991]  to  obtain  the 
mean  square  scattered  pressure  for  a  discrete  random  medium.  A  correlation 
function  or  a  spectral  function  characterises  the  random  medium.  This  is  obtamed 
from  a  large  collection  of  temperature  data  and  will  be  a  function  of  the 
envirorunent.  The  covariance  N2  is  given  by 


N2(>-r')  =  li^N{r-r'). 


I 


where  the  correlation  coefficient  N  is  defined  as 


^  be  be 
,  (r,v)  -(r',v) 

N(r-r')=  - 

c 


cind  |J.  as 


f(r.v)] 


The  symbol  (  )  indicates  averaging  over  different  temperature  profiles,  V .  is 
the  mean  square  fractional  speed  fluctuation  at  position  r .  The  spectral  function 
is  given  by  where  F(K)  is  the  Fourier  transform  of  N{r  -  r'). 

K  =  {K^,K^,K, ),  where  K  depends  on  the  incident  wave  number  and  the 
direction  of  incident  and  scattered  beam.  In  spherical  coordinates, 

Kx  =fco(sin0,  cos<t»i  -sinSj  cos<t»s ) 

Ky  =  <:o(sin0,  sin((», -sin0s  sinctij )  (7) 

=  ^0  (cos0,  -COS(>5  ) 


where  /  is  the  direction  of  the  incident  beam,  and  J  the  direction  of  the  scattered 
beam.  0  and  (|>  are  the  polar  and  azimuthal  angles  respectively.  For  back- 


scattered  pressure,  K  =  2kf^. 

Two  forms  of  correlation  and  spectral  functions  have  been  used  to  deal  with 
temperature  fluctuations  that  occur  at  the  fine  structure  and  micro  structure  level. 
One  is  isotropic  that  is  characteristic  of  micro  structure  fluctuations,  and  the  other 
is  anisotropic  and  characteristic  of  fine  structure  fluctuations.  The  isotropic 


functions  F{K)  and  N{x,y,2)  are  given  by  Liebermann  [1951],  Chernov  [1960, 
plO], 


F(/C)  =  (- 


8(11) 


% 


)e 


N{x,y,z)  =  e 


(7a) 
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where  o;,  is  the  correiadon  length. 

The  anisotropic  functions  F{K)  and  N(x,y,z)  are  given  by  Unni  and  Kaufnaan 
[19811, 


F(K)=[ 


P^20 


-\zU  -(x’^yV 

N(x,y,z)  =  e 


(7b) 


where  Zg  is  the  correlation  distance  in  the  depth  direction  and  is  the 
correlation  distaiKe  in  the  radial  direction.  The  form  of  the  anisotropic  function  is 
obtained  from  a  statistical  model  where  the  thermal  fine  airucture  in  the  .  ertical  is 
a  steppy  Poisson  process  and  where  the  layered  fine  structure  is  advected 
horizontally  and  vertically  by  internal  waves.  The  basis  of  the  steppy  Poisson 
process  is  the  experimental  discrete  probability  function  of  the  observed  layer  size 
[Joyce  and  Desaubies,  1977|.  The  introduction  of  a  steppy  Poisson  process  in  the 
vertical  introduces  the  exponential  function  in  the  anisotropic  correlation  function. 
However  it  must  be  noted  that  the  exponential  function  introduces  a  peculiarity, 
namely  its  partial  derivative  with  respect  to  Z  at  the  origin  differs  from  zero.  This 
implies  that  die  refractive  index  fluctuation  is  a  discontinuous  function  [Chernov, 
19^,  p9J.  The  exponential  function  is  commonly  used  however  because  it  often 

yields  realistic  results.  The  parameters  a^.z^.p^,  and  \l^  are  dependent  on  the 
environment. 

For  a  medium  where  the  mean  speed  of  sound  is  c,y  and  where  the  sound 
source  is  continuous  and  monochromatic,  the  mean  square  scattered  pressure  at  a 
distance  r ,  very  much  greater  than  the  dimensions  of  the  scattering  region  of 
I  olume  V,  is  given  by  (Chernov,  1960.  p51) 


|Ps(f)l^ 


F(k) 


(8) 


where  Icq  =  — /  W  ,  are  the  angular  frequency  and  amplitude  of  the  sound 
Co 

wave  respectively. 

Assuming  the  initial  direction  of  the  incident  beam  to  be  along  the  X  axis, 

0,  =  90*  and  <|)j  =  0’  let  us  evaluate  the  term 


F,(K)  = 


for  the  anisotropic  function  and 


(9) 


(10) 
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for  tfie  isotropic  function.  The  function  /[(/T)  contains  all  the  angular  and  the 
frequency  dependence  of  ttie  mean  square  scattered  pressure.  Fig.  7(a)  is  a  plot  of 
the  anisotropic  I\  (eqns  7,  9)  as  a  function  of  0^  at  1  MHz  for  two  scattering 
directions  tt,  =  89.8“  and  89.6'.  /\  is  plotted  on  a  logarithmic  scale  with 
=  74  m,  and  z^=2  m  (Unni  and  Kaufman.  1981  and  Joyce  and  Desaubies, 

1977],  These  lengths  were  obtained  for  the  Sargasso  sea,  and  are  not  universal 
parameters.  From  F^.  7(a)  we  observe  that  in  going  from  89. 8°  to  89. 6“  there  is  a 
lai^^e  decrease  in  the  magrtitude  of  .  Also,  the  function  decreases  very  rapidly 
with  increasing  azimuthal  scattering  angle  0^ .  This  indicates  that  scattering  is 
confined  to  a  small  angle  close  to  the  direction  of  propagation. 

Fig.  7(b)  is  a  plot  of  the  anisotropic  (eqns  7,9)  for  a  particular  scattered 

direction  0^  =  89. 8“  as  a  function  of  0^,  for  the  two  different  frequencies  1  MHz 
and  100  kHz.  The  results  indicate  that  as  the  frequency  decreases  scattering 
spreads  over  a  wide  angle.  However  the  spread  is  still  small,  the  function 
decreasing  rapidly  with  irKreasing  0^ . 

Fig.  7(c)  is  a  plot  of  the  isotropic  F^  (eqns  7,10)  as  a  function  of  the  azimuthal 
scattering  angle  0^,  for  the  two  angles  0^  =  89.9'  and  89.8°,  with  =  0.  7  m. 
Changing  the  scattered  direction  from  the  incident  direction  of  propagation,  there 
is  a  marked  decrease  in  F^ .  However  the  width  of  the  scattered  beam  using  the 
isotropic  spectral  function  is  much  greater  than  for  the  anisotropic  function.  The 
above  results  iiKlicate  that  high  frequency  (X  «  a^)  scattering  is  strongly 
directional,  the  direction  being  close  to  the  direction  of  propagation,  JChemov, 
1960,  p52]. 

Next  the  evaluation  of  the  total  attenuation  both  for  the  isotropic  and 
anisotropic  functions  is  undertaken.  Starting  from  equations  (.SO),  (SI)  and  (52)  of 
Chernov  [1960,  pp.S3,  .541,  the  attenuation  (a)  is  given  by, 

K  2r  _  2 

a  =  JsinWeJd(p)p,(r)p-^.  (11) 

0  0  'So  V 


The  expression  for  |ps  (r)|^  is  given  by  eqn  (8).  For  the  isotropic  correlation 
fuTKtion  it  is  independent  of  ttte  direction  of  propagation  and  can  be  easily 
evaluated,  (Chernov,  1960,  p55]; 

a{dB  /  m)  =  4. 34ii^iko^flo  Vn  (1-  )  (12) 

At  high  frequencies  it  varies  as  k/.  This  arises  since  for  large  »  1,  and 

thus  e”"®  *®  -» 0 .  Fig.  8(a)  shows  the  attenuation  in  dB/m  (eqn  12)  at  the  two _ 

frequence  1  MHz  and  100  kHz,  as  a  function  of  for  =  5.1X10”*®.  This 

corresponds  to  a  root  mean  square  temperature  fluctuation  of  7 X10~^  C  deg. 
Although  widi  increase  in  frequency,  scattering  becomes  confined  to  a  small 
region  dose  to  the  direction  of  propagation,  the  increase  in  attenuation  at  small 
angles,  causes  the  total  attenuation  to  increase. 

In  dte  anisotropic  case,  attenuation  will  depend  on  die  direction  of  propagation. 
Where  the  direction  of  propagation  is  along  the  Z  axis,  the  attenuation  is  given  by 

— -  5  -0.25pg*Jsin’e 

a(dB/m)  =  4.34p*i:*p?Zo  I ‘iesine— s-s - - -  (13) 

i  lko*zj(l-cose)2+lj 
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XfifB) 

Figures:  (m)  Attenuatim  (o^  in  dB/m  (eqn.  12)  mt  1  )AHz  and  100  kHz  as  a  function  of 
%  >1^  =  5.1  X 10"^°.  (b)  AUemtation  (o)  in  dB/m  (ajn/43)  at  1  MHz  and  100  kHz  as  a 
function  t^ZQfarthe  anisotropic  function.  po  =  50m,  p.^\s.l  x  10’“ 
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Integration  over  angles  in  the  above  integral  arises  due  to  scattering  being  in  all 
directions.  The  integral  is  evaluated  numerically  using  the  NAG  routine  DOIGAF. 
Since  high  frequency  scattering  is  confined  close  to  the  direction  of  propagation, 
in  this  sector  a  fine  interval  has  to  be  employed  to  obtain  accurate  results  for 
attenuation.  To  obtain  a  relative  error  of  less  than  1%,  for  0  upto  0.1  radians, 
integration  intervals  were  found  to  be  of  10  ^  radians  and  10^*  radians  at 
frequencies  of  100  kHz  and  1  MHz  respectively.  Beyond  0.1  radians  a  much 
coarser  interval  is  used.  Fig.  8(b)  shows  the  attenuation  (eqn  13)  at  1  MHz  and  1(X) 
kHz  as  a  fuiKtion  of  the  vertical  correlation  length  for  =  50  m  and 

=  5. 1X10“*° .  The  magnitudes  are  of  the  same  order  as  the  isotropic  function 
and  show  a  dependence.  There  is  negligible  variation  with  p^.  With 
propagation  along  the  vertical  Z  direction,  and  attenuation  being  confined  to 
angles  close  to  the  direction  of  propagation,  the  radial  correlation  length  is 
expected  to  have  a  negligible  effect. 


4.?  Amplitude  and  Phase  Fluctuations  in  Wave  and  Ray 
Theory 

In  the  last  section  we  found  high  frequency  acoustic  scattering  confined  to  regions 
close  to  the  direction  of  propagation.  This  effect  indicates  amplitude  and  phase 
fluctuations  will  become  important.  This  section  looks  at  obtaining  estimates  of 
these  fluctuations.  Consider  an  experiment  where  the  receiver  is  in  the  direction 
of  propagation  and  separated  by  a  distance  L  from  the  source,  which  transmits  a 
continuous  wave  monochromatic  signal.  The  receiver  receives  the  direct  wave  as 
well  as  waves  scattered  by  the  media.  The  superposition  of  the  direct  and 
scattered  waves  gives  rise  to  fluctuations  in  the  received  signal.  The  Rytov 
method,  (Chernov,  1960,  p61]  is  used  to  obtain  estimates  of  the  mean  square  log 

amplitude  fluctuations  3nd  the  mean  square  phase  fluctuations, 

(A(t»)‘ .  Here  A.  ^4^  are  the  amplitudes  and  <*),  <j>Q  the  phase  of  the  received  and 
transmitted  signals  respectively.  A<()  equals  0  — .  The  Rytov  method  is 
applicable  when  the  relative  amplitude  change,  and  pha.se  change  over  a 
wavelength  is  small.  There  is  no  restriction  on  the  total  change  of  these  quantities. 
Consider  the  source  being  placed  at  (0, 0, 0)  and  the  receiver  at  (L,  0,  0).  Then, 
following  Chernov  (1960,  p70),  we  have: 


(ln(^))2=0.5koV^Ii-J2) 

(14) 

(A0)^=O.5k^ti^(Ii  +  l2), 

(15) 

where  /,,  /j  are  integrals  that  depend  on  the  correlation  coefficient  and  are 
defined  by  Chernov  (1960,  p70].  Let  us  evaluate  the  integrals  for  the  two 
different  correlation  coefficients. 

For  the  isotropic  correlation  function, 

1  (l-T) 

Jl  =  L*Jdx  J  dy 

0  -X 


(16) 


1  l-» 


l2  =  l}jdx  j  dy 


0  ~x 


.4L^(2-2x-y)y 

/»54 


(17) 


These  are  obtained  by  starting  from  equations  (118)  and  (119)  of  Chernov  [1960, 
pTO]  and  without  making  any  further  approximations  (such  as  L  »  or 
»  1). 

Making  the  approximation  that  L  »  and  »  1,  Chernov  obtained,  sec 
his  eqns  (159),  (160)  in  Chernov  [1960,  p83). 


(ln(— ))^  =0.5Viitp^*o«oL(l-  — arctan(D)) 
Ao  D 

(A0)^  =0.5-</KVL^kQaQL(\  +  ^iTct!n\{D)) 


(18) 

(19) 


where  D  = 


4L 

Mo 


In  the  limit  —4  «»  one  can  obtain  simple  analytical  results.  In  this  case, 

J, er/(— )  +  (e  /“o-i)^) 

vJt 


(20) 


where  erf{x)  is  the  error  function  of  argument  X .  For  »  1 ,  using  eqns  (16) 
and  (17), 


1  i-»  Cl/' 


-1  -  ^2  “  f  f  -  2Ar  -  1/)^  -  1 

*=5-0  i  % 


(21) 


In  the  limit  ^  /j  Oand  therefore. 


(A<»^ «•/(— )  +  (r  /  "O-D-^l 

Oq 


(22) 


One  can  also  obtain  this  result  from  ray  statistics,  using  the  expression  for  the 
mean  square  transit  time  fluctuation  given  in  Chernov  [1960,  p291. 
Substituting  the  expression  for  (/,  -I2)  (eqn(21)),  into  eqn  (14)  the  relative 
intensity  fluctuatkm  at  large  becomes. 


(ln(^))^-»^V^Jd*  Jdy  r  '/‘^[(2-2x-y)2-y2l  (23) 

0  -I 

There  is  no  restriction  on  the  value  of  the  horizontal  range  in  the  above  equation. 
When^  » 1,  eqn  (22)  becomes. 
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->  li^kQLjnaQ 


(24) 


Chernov  obtained  the  same  result  in  the  ray  approximation,  using  the  isotropic 
gaussian  correlation  function  and  evaluating  the  integral  given  in  Chernov  [1960, 
p29].  When  » 1,  eqn  (23)  reduces  to 


Eqn  (25)  is  the  same  as  that  obtained  by  Chernov  in  the  ray  approximation, 
[Qiemov,  1960,  p34].  Clearly  his  expressions  in  the  ray  model  are  therefore  valid 
otJy  when  L»aQ. 


If 


« 1,  eqn  (22)  reduces  to. 


(26) 

For  the  anisotropic  correlation  function,  equations  (118)  and  (119)  of  Chernov 
[1960,  p  70]  are  again  the  starting  point  for  the  evaluation  of  the  integrals  /j  and 
/j.  When  evaluated  without  invoking  any  simplifications,  one  obtains  for  the 
integral 

1  -Ly/ 

Ii=L^jdx  jdij  e  /fo[M(f>i,po)/C(6,,Zo)  +  /fi((^,,Po)/D(fc,,2o)]  (27) 

0  .V 


The  expression  for  the  integral  /,  is  given  by 


1  1-* 

l2=L^jdx  jdij  e  /Po[M(&2,Po)/C(fc2,2o)  +  fB(^’2«Po)^C)(f72,Zo)] 

0  -I 

where  b2=^L(2-2x-y)- 


(28) 


M,  IB,  IC  and  ID  are  given  by. 
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/C(b,Zo)  =  (cos(<j)(0.5-S(«j))-sin(q)(0.5-C(^/)) 

ID(b,Zo)  =  (cos(<j)(0.5-C(q))+sin(<,)(0.5-S((j))  (30) 


.  .  ,  0.25 

in  which  q  =  — s— . 

^  z^b 

S(q)  and  C(^)  are  the  Fresnel  sine  and  cosine  integrals,  where 

1 

0 

<) 

C(q)  =  J  cos(t)dt 

In  the  limit  — >  <»,/j  -/j  — »  0,  with  /,  being  the  same  as  eqn  (20)  except  that 

Pq  replaces  .  Chernov  has  not  treated  the  anisotropic  case.  Equations  (16),  (17), 
(27)  and  (28)  do  not  have  any  restriction  on  the  range  or  frequency,  and  are 
applicable  at  all  ranges  and  frequencies. 


4.3  Numerical  Estimates  of  Intensity  and  Phase  Fluctuations 
for  the  Isotropic  Correlation  Function  in  the  YJave  and  Ray 
Models 

In  this  section  the  isotropic  correlation  function  is  used  to  obtain  some  numerical 
estimates  for  intensity  and  phase  fluctuations  at  high  frequencies.  The  dependence 

of  the  estimates  on  the  isotropic  correlation  length  ,  frequency  (/),  and 

range(Z,)  is  examined.  Also,  wave  model  results  are  compared  with  the  results 
from  the  Chernov  model  and  the  statistical  ray  model.  A  brief  discussion  on  the 
trends  observed  in  the  simulation  results  is  given.  A  frequency  of  100  kHz  and  a 
range  close  to  100  m  are  often  used  in  target  strength  measurement  of  mines.  This 
is  the  reason  for  the  choice  of  these  values  for  frequency  and  range  in  most  of  the 
calculations  presented  in  this  section. 

In  Fig.  9(a),  comparisons  of  the  rms  relative  intensity  fluctuation  (eqns  (14, 16 
and  17))  and  the  rms  phase  fluctuations  (eqns  (15, 16  and  17))  are  made  with  the 
Chernov  result  (eqns  (18)  and  (19)),  as  a  function  of  range.  The  NAG  routines 
DOl  AJF  and  DOlGAF  are  used  for  the  numerical  evaluation  of  the  iimer  and  outer 
integrals  in  /,  and  f  respectively.  Using  a  relative  error  of  10'^  for  the  inner 

integral  and  an  interval  of  for  the  outer  integral,  a  relative  error  of  less 

than  1%  is  achieved  for  the  double  integral.  The  parameters  used  in  the 

compariscm  are  /  =  l(X)kHz,  O',,  =2  m  and  (i*  =5.1X10"^®.  The  value  of 

corresponds  to  a  rms  temperature  fluctuation  of  7X10~^C  deg  and 
kfyOQ  *>  838  » 1.  The  plot  shows  dut  for  L  greater  than  50  m,  ^ere  is  an  error  of 
less  than  5%  in  the  intensity  between  the  results  of  Chernov  and  the  wave  model . 
At  large  range  where  L  is  very  much  greater  than  both  models  agree.  In  Fig. 
9(b)  comparison  is  made  between  the  models  as  a  function  of  This  is  carried 

out  with  L  — 100  m,  /  =  100 kHz,  and  =  5.1X10"*®.  The  maximum  value  of 
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Oj  used  is  0. 75  m,  thus  L  »  a^.  The  plot  shows  that  the  two  models  agree  only 
for  »  80. 

Figure  10  compare  the  rms  relative  intensity  and  rtns  phase  obtained  from  the 
ray  model  (eqns  (23, 22))  with  that  of  the  Chernov  ray  m^el  (eqns(25, 24)),  as  a 

functionof  range  L.  The  parameters  used  are  /  =  200 kHz,  =  5.1X10“*®, 

^  =  2  m,  and  =  1676.  Figure  10(a)  compares  the  rms  relative  intensities 
and  Fig.  10(b)  compares  the  rms  phases.  The  graphs  indicate  that  the  Chernov  ray 
model  results  agree  widt  our  ray  model  results  with  an  error  of  less  than  5%  only 

for  L  greater  than  20  m.  At  this  value  of  L,  =10.  In  oUrer  words  the 

Chernov  ray  model  (eqns  (25, 24)),  is  valid  only  for  L  very  much  greater  than  . 

In  Fig  11(a)  the  rms  relative  intensity  fluctuations  (eqns  (14, 16, 17)),  and  in  Fig. 
11(b)  the  rms  phase  fluctuations  (eqns  (15, 16, 17)),  are  plotted  for  different 
isotropic  correlation  lengths  Oq,  with  /  =  100  kHz,  L  =  100  m,  and 

=  5. 1X10“*® .The  correlation  length  varies  from  10  cm  to  2  m.  The  integrals 
(eqns  (16, 17))  are  evaluated  numerically  as  before.  The  maximum  variation  in  the 
rms  relative  intensity  for  this  range  of  is  less  than  0.3  dB  and  the  maximum 
phase  fluctuation  about  0.15  radians.  The  peak  observed  around  0.5  m  in  the  rms 
intensity  curve  is  interesting,  and  this  is  discussed  later,  along  with  the  trend 
observed  in  the  fluctuations  as  a  function  of  . 

Figure  12(a)  shows  the  variation  of  the  rms  relative  intensity  fluctuations  (eqns 
(14, 16, 17))  for  different  .  The  variation  in  corresponds  to  a  variation  in  the 

rms  temperature  fluctuation  from  7X10“^  to  3X10“^  C  deg  at  a  particular  depth 
and  salinity.  Fig.  12(b)  shows  the  variation  in  the  rms  phase  (eqns  (15, 16, 17)).  The 
rms  relative  intensity  fluctuation  varies  by  about  1.2  dB  and  the  rms  phase 
fluctuation  by  about  0.3  radians  over  titis  interval  of  temperature  fluctuation.  The 
calculations  are  carried  out  with  /  =  l(X)  kHz,  =  0. 7  m  and  Z,  =  100  m.  We 
have  plotted  the  results  here  on  a  semi-logarithmic  scale  due  to  the  range 

involved  in  .  However  if  we  plot  it  on  a  linear  scale  against  (n*  )® a  linear 
dependence  is  observed  with  the  rms  relative  intensity  and  phase  fluctuations,  in 
accordance  with  eqns  (14)  and  (15)  respectively. 

Figures  13(a)  and  (b)  show  the  variation  of  die  rms  relative  intensity  fluctuations 
(eqns  (14, 16, 17))  and  rms  phase  fluctuations,  (eqns  (15, 16, 17))  as  a  function  of 
frequency,  where  the  frequency  varies  from  100  kHz  to  400  kHz.  The  parameters 

used  in  the  calculation  are  =  5.1X10“*®,  =  0. 7  m  and  L  =  100  m.  The  rms 

relative  intensity  fluctuations  initially  show  a  non  linear  variation  with  frequency, 
but  as  the  frequency  increases,  a  linear  dependence  is  seen.  Simulations  carried 
out  at  high  frequencies  show  that  the  relative  intensity  approaches  the  result  of 
eqn  (25),  where  it  will  be  independent  of  frequency.  This  is  expected  sirKe 
I. »  dg  in  this  example.  The  rms  phase  fluctuation  on  the  other  hand,  exhibits  a 
linear  variation  even  in  the  frequency  range  of  100  kHz  to  400  kHz  (it  follows  the 
trend  predicted  by  either  eqn  (22)  or  (24)).  This  indicates  that  for  phase 
fluctuations,  agreement  with  the  ray  result  can  be  obtained  at  a  lower  frequency. 

Figures  14(a)  and  (b)  show  the  variatkm  of  d\e  rms  relative  intensity  fluctuations 
(eqns  (14, 16, 17))  and  rms  phase  fluchiations  (eqns  (15, 16, 17))  as  a  function  of 
range  L,  where  the  range  varies  from  50  m  to  300  m.  Odier  parameters  taken  for 

die  calculations  are  |i*  =  5. 1X10“*®,  =  0. 7  m  and  /  =  100  kHz.  The  rms 

relative  intensity  shows  a  near  linear  relationsh^).  In  the  ray  approximation  it 
varies  as  I)  (eqn  (25)),  showii^  that  die  frequency  has  to  be  much  higho'  for  this 
limit  to  be  reacl^.  The  rms  phase  fluctuations  initially  show  a  non  linear 


behaviour  of  becoming  linear  with  increase  in  range.  In  the  ray 
approximation  the  variation  with  rai^e  is  non  linear  and  is  of  the  form  (eqn 
(24)).  It  will  be  seen  later  that  for  a  particular  hequency,  the  smaller  the  range,  the 
better  is  the  agreement  with  the  ray  result. 

Figure  15(a)  compares  the  rms  phase  fluctuations  obtained  from  the  stochastic 
ray  model  (eqn  (22)),  with  the  wave  model  (eqns  (15, 16, 17)),  as  a  function  of 

^QdQ.The  calculations  are  carried  out  with  L  =  100  m  and  =  5.1X10“^*^.  The 
differences  between  the  two  models  decrease  with  increase  in  k^a^.  Since  the 
stochastic  ray  model  is  obtained  from  the  wave  model  in  the  limit  the 

observed  results  follow  the  expected  trend. 

The  correlation  distance  is  a  measure  of  the  decay  length  of  the  coherence  of 
fluctuations.  The  bigger  the  length  the  larger  the  d^tance  at  which  fluctuations 
become  uncorrelated.  For  a  range  greater  than  a^,  the  larger  the  value  of  Oy,  more 
time  is  spent  in  going  through  the  coherent  region  of  fluctuations.  This  increases 
the  phase  fluctuations,  as  seen  also  from  eqn  (24).  On  the  other  hand  if  the  range  is 
mu^  smaller  than  we  would  expect  the  phase  fluctuations  to  be  independent 
of  Uq  ,  as  seen  also  from  eqn  (26). 

Variation  of  mtel^sity  fluctuation  with  correlation  length  is  more  involved  than 
the  variation  of  phase  fluctuations.  To  understand  this  let  us  consider  Fig.  1.  The 
diffraction  or  the  size  parameter  A  is  given  by  [Flatte,  1979,  eqn.6.2.3J, 


A  = 


L 

bi^kQ 


(31) 


The  strength  parameter  <1>  is  the  rms  phase  fluctuation  in  the  ray  approxintabon 
and  is  given  by  eqn  (22).The  geometrical  approximation  or  the  ray  model  is  valid 
when  diffraction  effects  are  small,  that  is  A  <  1  and  A4>^  <  1 .  The  Rytov  method 
is  applicable  for  A  >  1  but  AO'  <  I .  Fig.  16(a)  is  a  plot  of  the  diffraction 
parameter  (eqn  (31)),  as  a  function  of  a^,  for  f  =  100  kHz  and  Z,  =  100  m.  The 
plots  indicate  that  diffraction  becomes  Lnxportant  for  less  than  0.5  m.  We 
observed  before  that  the  maximum  in  the  rms  intensity  as  a  function  of 
occurred  close  to  0.5  m.  The  Fresnel  zone  radius  is  defined  as  the  size  of  the 
inhomogeneities  that  will  cause  the  scattered  path  length  from  source  to  receiver 
to  differ  from  the  unscattered  path  length  by  half  a  wavelength.  When  this 
happens  there  will  be  interference  between  the  direct  and  scattered  rays.  At  100 
kHz  over  a  path  length  of  100  m,  this  is  evaluated  and  is  equal  to  0.5  m.  This  is 
close  to  the  value  of  at  which  the  maximtmt  in  the  rms  relative  intensity 
occurred  (see  Fig.ll(a)). 

Figure  16(b)  is  a  plot  of  O  the  rms  phase  fluctuation  obtained  in  the  ray 
approximation,  given  by  eqn  (22),  as  a  function  of  .  Figxue  16(c)  is  a  plot  of 

A^^  as  a  fuiKtion  of  Oq  .  Correlation  length  {a^  )  varies  from  10  cm  to  2  m.  Other 

parameters  used  in  the  calculations  are  =  5.1X10~^*^,  /=  100  kHz  and 

L  =  100  m.  Tfre  results  indicate  that  we  are  in  the  unsaturated  region  of  the 
I  A  -  4>  space.  Unlike  the  nns  phase  fluctuations,  intensity  fluctuations  decrease 

J  with  correlation  distance.  SiiKe  diffraction  effects  decrease  with  increase  in 

correlation  leiigth,  the  ray  model  becomes  a  good  approximation  to  evaluate 
j  fluctuations  at  large  correlation  length.  In  the  geometric  approximation  intensity 

i  fluctuations  decrease  with  increase  in  the  conation  length  [Chernov,  1960,  p34]. 

^  We  would  therefore  expect  the  same  trend  in  the  results  of  intensity  fluctuations 

from  the  wave  model  litis  trend  is  seen  in  Fig.ll(a). 
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Figure  10:  Comparison  of  the  (a)  rms  relative  intensity  and  (b)ms  phase  obtained  from  the  ray  model,  eqns.  (23,  22)  unth  that  of  the  Chernov  ray 
model,  eqns.  (25,24)  as  a  function  of  the  range  L.  f=  200  kHz,  =  5.1  x  10"^°,  aQ  =  2m,kao  =  1676. 
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(14, 15, 27-30)  for  the  anisotropic  medium. 


and  eqns.  (14, 15, 27-30)  for  the  anisotropic  medium. 
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Figure  13:  RMS  relative  intensity  (dB)  and  rms  phase  (radians)  as  a  function  of  frequency  for  the  isotropic  (a,  b)  and  anisotropic  (c,  d)  media. 
5.1  X  L  =  100  m.  For  the  isotropic  media,  Aq  =  10  m,  and  for  the  anisotropic  media  zq  =  1  m,  pg  =  50  m.  Used  eqns.  (14-17),  for  the 

isotropic  medium  and  eqns.  (14, 15,  27-30)  for  the  anisotropic  medium. 


'.ity  (dB)  and  rms  phase  (radians)  as  a  function  of  the  range  (L)  for  the  isotropic  (a,  b)  and  anisotropic  (c,  d)  media. 

For  the  isotropic  media,  ag  =  I.O  m,  and  for  the  anisotropic  media  Zg  =  I  m,  pg  =  50  m.  Used  eqns.  (14-17),  for  the 
15,  27-30)  for  the  anisotropic  medium. 
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♦  Figure  16:  (a)  Plot  of  A  (eqn  (31)),  (b)  Plot  of<P  (eqn  (22)),  (c)  Plot  of  A<D^  as  a 

function  of  the  isotropic  correlation  length  Oq.  f=  100  kHz,  L  =  100  m,  =  5. 1  x  JO*®. 
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4.4  Numerical  Estimates  of  Intensity  and  Phase  Fluctuations 
for  the  Anisotropic  Correlation  Function  in  the  Wave  and  Ray 
Models 

This  section  gives  some  numerical  estimates  of  high  frequency  intensity  and 
phase  fluctuatior«  using  the  anisotropic  correlation  function.  The  variation  with 

the  radial  and  vertical  correlation  lengths,  range,  p' ,  and  frequency  are 
examined.  In  addition,  the  wave  model  results  are  compared  with  the  results  of 
the  stochastic  ray  model. 

Figures  11(c)  and  Fig.  11(d)  show  the  rms  relative  intensity  fluctuations  (eqns 
(14, 27-30)),  and  rms  phase  fluctuations  (eqns  (15, 27-30)),  for  different  values  of 
the  vertical  correlation  length  .  As  for  the  isotropic  case,  the  NAG  routines 
EWl AJF  and  DOlGAF  are  used  to  evaluates  the  integrals  /,  and  /,  (eqns  (27)  and 

(28)) .  Using  an  interval  of  5X 10”^  for  the  outer  integral,  a  relative  error  of  less 
than  1%  is  achieved  in  the  double  integral.  The  NAG  routines  S20ADF  and 
S20ACF  evaluate  the  Fresnel  integrals  S{q)  and  C{(j)  (eqn  (30))  respectively  .  The 
correlation  length  varies  from  10  cm  to  2  m.  The  values  of  the  other  parameters 

used  are  L  =  100  m,  =  50  m,  =  5. 1X10“^®,  and  /  =  100  kHz.  The  results 
using  the  anisotropic  correlation  function  are  in  general  similar  to  that  obtained 
using  the  isotropic  correlation  function,  except  that  the  variations  are  much 
higher.  The  variation  in  the  relative  intensity  is  around  1.25  dB  compared  to  about 
0.3  dB  for  the  isotropic  case.  For  the  phase  fluctuations  the  variations  are  about  0.4 
radians  as  compared  to  0.15  radians  for  the  isotropic  case.  This  increase  could  be 
due  to  increased  scattering  at  angles  close  to  the  direction  of  propagation,  as  a 
result  of  the  anisotropy. 

Figures  17(a)  and  (b)  show  the  rms  relative  intensity  (eqns(14,  27-30)),  and  phase 
fluctuations  (eqns  (15,  27-30)),  for  different  values  of  the  radial  correlation 
distance  p^.  The  radial  correlation  distance  varies  from  10  cm  to  100  m  The 

values  of  the  other  parameters  used  are  Zq  =  1  m,  =  5. IX 10"“^,/  =  100  kHz, 

«md  L  =  100  m.  The  relative  intensity  variations  and  phase  variations  are 
approximately  1.5  dB  and  0.8  radians  respectively. 

Figures  12(c)  and  (d)  show  the  rms  relative  intensity  and  phase  fluctuations  for 

different  values  of  p,‘ .  The  variation  in  p'  corresponds  to  a  temperature 

fluctuation  that  varies  from  7  X  10~^  to  3X10”^  C  degree  at  a  particular  depth  and 
salinity.  The  vedues  of  the  other  parameters  used  are  /  =  1 00  kHz,  L  =  100  m, 

Zq  =  0.7  m,  and  pg  =  50  m.  The  variation  in  acoustic  fluctuation  is  much  higher 
than  that  obtained  with  the  isotropic  correlation  function.  The  variation  in  the  rms 
relative  intensity  is  about  6  dB,  and  the  variation  in  phase  is  about  2.4  radians, 

over  the  range  of  p*  used.  This  shows  that  the  rms  intensity  and  phase 
fluctuations  are  sensitive  to  p' .  The  plot  shown  is  a  semi-logarithmic  plot  due  to 
the  range  involved  in  p' .  However  a  plot  of  the  acoustic  fluctuations  against 
(p^  ’  on  a  linear  scale  gives  a  straight  line  that  passes  through  the  origin. 
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Figure  17;  (a)  RMS  relative  intensity,  eqns  (14,  27-30)  (b)  and  rms  phase  eqns  (15, 27- 
30)  as  a  function  of  the  radial  correlation  length  po  of  the  anisotropic  medium, 

f  =  100  kHz,  L  =  100  m,  =  5.1  x  lO"*®  and  Zq=  1  m. 
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Figures  13(c)  and  (d)  show  the  variation  of  rms  relative  intensity  and  phase 
fluctuations  as  a  function  of  frequency,  where  the  frequency  varies  from  100  kHz 

to  400  kHz.  The  values  of  the  other  parameters  used  are  =  5. 1 X 10”'° ,  L  =  100 
m,  Pq  =  50  m  and  =0.7  m.  The  rms  relative  intensity  shows  a  near  linear 
variation.  The  rms  phase  fluctuation  increases  linearly  with  frequency,  as 
predicted  by  eqn  (22).  The  result  is  similar  to  that  obtained  for  the  isotropic  case. 
Agreement  widi  the  ray  result  is  obtained  at  a  lower  frequency  for  the  rms  phase, 
than  for  the  relative  intensity.  Comparison  of  the  results  for  the  two  media 
indicates  that  for  the  anisotropic  medium  the  rms  phase  is  about  6  times  larger 
and  the  relative  intensity  about  5-12  times  larger  than  the  isotropic  medium,  for 
the  parameters  chosen  in  this  example. 

Figures  14(c)  and(d)  show  the  variation  of  rms  relative  intensity  and  phase 
fluctuations  as  functions  of  range,  where  the  range  varies  from  50  m  to  300  m.  The 

values  of  the  other  parameters  used  are  =  5. 1  X10“^° ,  Zq  =  0. 7  m,  =  50  m 
and  f  —  100  kHz.  The  rms  relative  intensity  shows  a  near  linear  relationship  with 
range,  with  values  on  the  average  approximately  5  times  larger  than  for  the 
isotropic  case.  The  rms  phase  variation  shows  a  non  linearity  of  Z,”  "  initially,  close 
to  that  predicted  by  the  ray  result  (eqn  (22, 24)).  As  the  range  increases,  the 
variation  becomes  linear  showing  that  for  a  particular  frequency,  the  shorter  the 
range,  closer  the  agreement  with  the  ray  result.  This  is  understood  on  the  basis  of 
the  parameters  A  and  d>'  introduced  for  the  isotropic  case,  both  of  which 
increase  with  range.  The  rms  phase  values  are  higher  by  5-7  times  the  value  of  the 
isotropic  case,  for  the  parameters  used  in  this  example. 

Figures  15(b)  and  (c)  compares  rms  phase  fluctuations  obtained  from  the  wave 
model  (eqns  (14, 27-30)),  with  that  of  the  stochastic  ray  model  (eqn  (22),  with  p^ 
replacing  a^).  This  is  shown  in  Fig.  15(b)  as  a  function  of  and  in  Fig.  15(c)  as 
a  function  of  ^oPo-  The  values  of  the  other  parameters  used  are 

p^  =  5.1X10"^‘’,L  =  100  m,and  /  =  100 kHz.  In  Fig.  15(b),  Po  =  50  m,  and  in  Fig. 
15(c),  Zq  =  1  m.  In  both  cases,  the  graphs  indicate  that  agreement  with  the  ray 
model  improves  with  increase  in  and  Ar^p^ .  Since  the  ray  model  results  are 
obtained  by  making  »  «>,  the  results  follow  the  expected  trends. 


5.  Phase  Structure  Function  using  Stochastic 

Ray  Theory 

In  the  previous  sections  we  obtained  estimates  of  the  variance  of  the  relative 
intensity  and  phase  fluctuations  at  a  point  in  space.  On  comparing  the  results  of 
phase  fluctuations  obtained  from  the  wave  and  ray  models,  it  can  be  seen  that  for 
L  =  100  m,  if  >  500,  the  difference  between  the  models  is  less  than  5%.  This 
can  be  seen  from  Fig.  15(a).  To  obtain  the  same  agreement  in  the  anisotropic 
medium,  IcqZq  and  liQPg  must  be  greater  than  600,  as  can  be  seen  from  Figs.  15(b) 
and  (c).  In  high  frequency  mine  imaging,  the  range  involved  is  a  few  metres  and 
frequencies  are  in  the  mega  hertz  band.  To  examine  the  accuracy  of  the  stochastic 
ray  model  to  predict  the  rms  phase  fluctuations  under  these  conditions,  wave 
model  calculations  (eqns  (15-17)  for  the  isotropic  medium  and  eqns  (15, 27-30)  for 
the  anisotropic  medium),  and  ray  model  calculations  (eqn  (22))  are  carried  out 

with  =  5.1X10“'°,/=  1  MHz ,  and  Z.  =  3  m.  The  results  of  these  calculations 
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are  shown  in  Fig.  18(a)  as  a  function  of  the  isotropic  correlation  length  a^,  and  in 
Fig.  18(b)  as  a  function  of  the  anisotropic  correlation  length  p^.  In  Fig.  18(b) 
Zq=2  m.  The  curves  in  Fig.  18  indicate  that  for  /  =  1  MHz,  and  L  =  3  m  the 
results  of  the  two  models  are  indistinguishable  graphically.  This  shows  that  for 
these  situations  the  stochastic  ray  model  can  provide  good  estimates  of  the 
variance  of  the  phase  fluctuations.  This  being  the  case,  one  may  proceed  to  obtain 
estimates  of  the  phase  structure  function  in  the  ray  approximation. 

The  phase  structure  function  (PSF)  is  defined  by  [Flatte  et  al  1979,  pl08] 


PSF(x)  =  l(HRi)-^(Ri))-(^(R2)-^(R2))f  (32) 


where  0(Ri ),<(>( R2)  are  the  phases  of  the  signal  at  receivers Ri,R2,  4)(R2),<|)(R2) 
their  mean  values,  and  x  the  distance  of  separation  between  the  receivers.  The 
phase  structure  function  gives  the  variance  of  the  phase  difference  between  two 
points  caused  by  sound  speed  fluctuation.  Expanding  PSF(x),  it  can  be  written  as 


PSF(x)  =  (A0(R,  )f  +  (  A0)(R2  ))2  -  2A(D(Ri  )A(|>(R2 )  (33) 


where  A(t>(R)  =  4»(R)-<1)(R).  lneqn(33)  (A(l>(Rj))^  is  the  variance  of  the  phase 
fluctuation  at  receiver  1  and  (A4i(R2  ))^  is  the  variance  of  the  phase  fluctuation  at 
receiver  2.  The  term  A(ti(R|  )A(t>(R2)  is  the  spatial  auto  correlation  function  of  the 
phase  fluctuations.  The  phase  structure  function  therefore  involves  the  variances 
of  the  phase  fluctuation  at  receivers  1  and  2  and  the  spatial  auto  correlation 
function  of  the  phase  fluctuations.  If  the  variance  of  the  phase  difference  between 
the  two  receivers  is  large  then  the  coherence  between  the  two  signals  is  small  and 
vice  versa,  {Esswein  and  Flatte,  1981).  Since  the  phase  structure  function  is  a 
measure  of  variance  of  the  phase  difference  at  the  two  receivers,  it  is  a  useful 
parameter  in  the  evaluation  of  coherence.  The  phase  structure  function  is  also  seen 
as  a  measure  of  correlation  between  the  signals  at  the  two  receivers.  If  the  signals 
are  completely  independent  then  PSF(x)  approaches  a  value  that  is  the  sum  of  the 
variance  of  the  phase  fluctuations  at  the  two  receivers. 

Let  us  now  estimate  the  phase  structure  function  due  to  thernnal  fluctuations  in 
the  ray  approximation.  Suppose  we  have  the  source  at  (0, 0, 0),  receiver  at 

(jR,  0, 0)  and  at  (A'j ,  0,  r ),  where  Xl  +  Z^  =  That  is  both  receivers  are  at 

the  same  distance  from  the  source  but  are  in  different  directions,  one  being  at  a 
different  depth.  Using  the  expression  for  PSF  in  terms  of  the  correlation  function, 
[Esswein  and  Flatte,  1981]  and  sustituting  the  isotropic  correlation  function  [eqn 
7(a)],  one  obtains 


where 


_  -rV 

PSf(z)  =  fcoV^>/5^"o[2Rer/(^^)+2^(e  Ao_i)_p„j  (34) 


(35) 


I 
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Figure  18:  Comparison  of  the  RMS  phase  fluctuations  from  the  wave  model,  for  an 
isotropic  (eqns  15-17)  and  anisotropic  (eqns  15,  27-30)  medium,  with  the  stochastic  ray 

model,  eqn  (22).  f  =1  MHz,  L  =  3m  and  =5.1  x  10"^®.  (a)  as  a  function  of  Oq,  the 
isotropic  correlation  length,  (b)  as  a  function  of  pQ,  the  anisotropic  radial  correlation 
length,  with  the  vertical  correlation  length  Zq  =  2  m.  The  2  curves  are  indistinguishable 
within  graphical  accuracy. 
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For  the  anisotropic  correlation  function  [eqn  (7b)], 

-iiV 

PSF{z)  =  k^V}^,U^Rerf(jy^)^^,{e  (36) 


where 

,  -Xl/  JhXW 

P22  =  (l+m2)(2(bpo  +  X2)  Je-<' ^d/  +  Po(<r  /P»  (37) 

b 


P33  =  R f  dx  r  Ier/(-5-^)  +  er/(^)]  (38) 

i  Po  Po 


in  which 


The  above  expressions  for  the  phase  structure 


function  are  obtained  assunung  the  path  of  the  ray  between  the  source  and 
receiver  to  be  a  straight  line.  This  is  a  reasonable  assumption  for  a  range  of  a  few 
metres  and  when  the  mean  gradient  is  small.  There  is  no  restriction  on  the 
separation  between  receivers. 

When  the  two  rays  are  almost  parallel  to  the  X  axis,  with  the  receiver  JR^ 
positioned  at  (R,  0, 2),  for  the  isotropic  correlation  function. 


PSf(2)  =  2«b?(l-e 


(39) 


and  for  the  anisotropic  correlation  function. 


PSf(2)  =  2«I>f(l 


-,%i 


(40) 


where  Of  is  the  mean  square  phase  fluctuation  at  ,  d>i  =  { A<t>(  Rj  ))^ ,  given  by 
eqn  (22).  For  the  anisotropic  function  p^,  replace  in  eqn  (22).  The  assumption 
that  two  rays  are  parallel  would  be  a  reasonable  approximation  at  large  range  and 
for  snull  separation  between  the  receivers.  One  can  also  derive  eqn(39)  in  the  ray 
approximation  for  an  isotropic  medium  using  R  »  Og ,  z  •>  og  >  Oe  spatial  auto 

correlation  coefficient  of  the  phase  fluctuations  obtained  by  Chernov  [1960,  pl07]. 
This  indicates  that  the  closed  form  eqns  (39)  and  (40)  are  valid  only  at  large  range 
and  for  small  separation  between  the  receivers.  No  expressions  are  given  by 
Chernov  for  small  R  or  for  die  anisotropic  case. 

Let  us  now  examine  the  results  of  some  numerical  simulations  carried  out  for 

the  PSF.  Consider  first  the  isotropic  medium.  Fig.  19(a)  is  a  plot  of 


PSF(z\ 


%\ 


(eqns  (34, 35, 22))  as  a  function  where  R  is  equal  to  3  m.  The  two  plots 


correspond  to  =1  and  2  m.  The  expression 


^^^(^^2 ’s  independent  of 


and  frequency,  and  depends  only  on  R,  and  z.  This  can  easily  be  seen  from 
eqns  (34, 35, 22).  Since  R  is  comparable  to  use  of  eqn  (39)  to  obtain 


would  be  inappropriate.  The  integral  diat  occurs  in  eqn  (35)  is  evaluated 
numerically  using  the  NAG  routine  DOl  AJF.  A  relative  error  of  less  dian  1%  in  die 
integral  is  ^iain^  by  specifying  the  relate  error  in  die  routine  to  be  less  dian 
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10~*.  A  qualitative  understanding  of  the  variation  with  is  obtained  using  eqn 
(39).  This  indicates  that  for  small  the  value  of  Z  to  reach  saturation  is  also 
small. 

Figurel9(b)  shows  the  plot  of  ®  function  of  where  Aq  =  0. 1  m 

and  R  =  3  m.  PSF(z)  is  calculated  using  eqns  (39, 22)  as  well  as  eqns  (34, 35).  For 
R  »  Oq,  eqn  (39)  would  be  appropriate  only  at  short  separations.  The  calculations 
are  carried  out  with  receiver  Ry  being  positioned  at  (R,0,z).  At  small 
separations,  where  z  <  O.Saq,  the  results  of  both  approximations  agree.  However, 
at  large  Z ,  discrepancies  appear,  where  the  parallel  approximation  ((eqn.(39)) 
fails.  At  very  large  vertical  separations,  both  approximations  give  a  value  of  2  for 

is  true  always  for  an  isotropic  medium  and  cannot  be  used  to 

distinguish  between  the  approximations.  For  an  isotropic  medium  the  first  two 
terms  of  the  PSF  given  in  eqn  (34)  are  equal  for  all  separations  z,  provided  the 
distances  between  the  source  and  receivers  are  the  same.  The  third  term  decreases 
to  zero  at  large  separatiorrs.  This  accounts  for  the  result  observed  at  large  vertical 
separations.  Figures  19(a),  (b)  show  that  as  decrease,  the  separation  at  which 
saturation  occurs  also  decrease.  As  decreases,  the  distance  over  which  the 
acoustic  phase  fluctuations  are  correlated  decreases,  which  in  turn  causes  the 
distance  at  which  saturation  occurs  to  decrease. 

Let  us  now  consider  an  anisotropic  medium.  In  this  case  the  phase  fluctuations 
at  positions  that  are  separated  from  the  source  by  the  same  distance  R  but  are  in 
different  directions  in  the  xz  plane  will  not  be  the  same.  When  the  angle  of 
inclination  of  the  receiver  measured  from  the  horizontal  direction  increases,  there 
is  a  decrease  in  the  variance  of  the  fluctuations  at  the  receiver.  This  is  because  the 
temperature  fluctuation  decays  much  faster  in  the  vertical  direction  than  in  the 
radial  direction,  resulting  in  a  reduction  in  the  acoustic  fluctuation.  For  large 

separation  we  would  therefore  expect  that  a  plot  of  a  function  of  ^ 

to  approach  1.  The  parallel  ray  approximation  (eqn. (40))  clearly  fails  at  large  Z 
where  it  approaches  a  value  of  2.  This  is  because  the  two  paths  to  the  receivers  are 

the  same  in  the  parallel  ray  approximation.  Figure  20  is  a  plot  of  ® 

/<Pi 

function  of  ^ .  The  calculations  are  carried  out  with  pg  =  75  m,  R  =  3 ,  and 
Zq  =  0. 1, 0. 5  m.  Since  R  «  pg,  eqns  (36,37,38,22)  are  used  to  evaluate  the  PSF.  The 
integral  that  appears  in  eqns  (37)  and  (38)  is  again  evaluated  numerically  using  the 
NAG  routine  DOIAJF.  A  relative  error  of  less  than  1%  for  the  integral  is  obtained 
by  specifying  the  relative  accuracy  to  be  less  than  lO”*  in  the  routine.  When  Zq  is 
small,  saturation  is  reached  at  small  vertical  separations.  The  results  display  the 
expected  trends. 


Figure  19:  Plots  of(PSF(z)/^)  as  a  function  of  (a)  z/R,  (b)  zla^for  an  isotropic  medium, 
z  is  the  vertical  separation  between  the  2  receivers,  and  ^  the  mean  square  phase 
fluctuation  at  receiver  1.  R  =  3m.  (a)  Curve  (1)  is  for  aq  =  2  m  and  curve  (2)  foraQ  = 

1  m.  Both  curxxs  plotted  using  eqns  (22. 34. 35).  (b)  Curve  (1)  plotted  using  eqns  (22, 

34, 35)  and  curve  (2)  using  eqn  (39).  In  (b)  both  curves  use  Oq  =  0.1  m.  Graph  (b) 
compares  the  linear  and  parallel  approximations. 
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Figure  20:  Plots  of  (PSF(z)/<t>i  )  (eqns  22,  36-38)  as  a  fiinction  ofz/R  for  an  anisotropic 
medium,  Zq  =  0.1  and  0.5  m.  K  =  3  m,  po  *  75  m. 


6.  Summary 

In  this  paper  suitable  models  to  obtain  high  frequency  acoustic  scattering 
estimates  due  to  temperature  changes  were  examined.  Previously  models  were 
restricted  to  long  range^  and  did  not  treat  the  anisotropy  of  the  fluctuations.  Iliere 
was  thus  a  need  to  correct  this  problem.  CXher  questions  that  were  looked  into 
include  the  validity  of  the  stochastic  ray  model  at  high  frequencies  (  say  100  kHz 
and  above),  theoretical  estimates  for  the  phase  structure  function,  and  estimates  of 
intensity  and  phase  changes  due  to  the  mean  temperature  gradient. 

Section  2  shows  some  experimental  temperature  and  temperature  gradient 
profiles.  They  show  the  type  of  data  obtained  experimentally,  and  give  a  feeling 
for  the  magrytudes  of  the  tenqxual  and  spatial  variation  of  temperature.  Near 
surface  layers,  there  are  temperature  inversions,  and  the  average  temporal 
varuLtkm  of  temperature  ranges  from  10'^  to  10'‘  C  deg  per  min  (Figs>(3, 4)).  Hie 
tenqseniture  gradient  is  respon^le  for  acoustic  refiraction.  The  average 
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temperature  gradient  in  the  vertical  direction,  obtained  from  the  ocean  data  near 
PerA,  shows  a  value  of  -0.1  C  deg  /m.  However  at  some  positions  in  the  profile 
the  temperature  gradients  are  much  higher,  close  to  -0.15  C  deg/m.  The  gamma 
probability  distribution  function  [Hayes  et  al,  1975]  fits  quite  well  to  the  discrete 
probability  function  ot  the  temperature  gradient. 

Section  3  gives  results  for  intensity  and  phase  changes  due  to  the  mean 
temperature  gradient.  Expressions  for  intensity  changes,  eqns  (1, 3),  and  phase 
changes,  eqns  (5, 6),  due  to  dte  mean  temperature  gradient  are  given.  They  were 
obtained  from  deterministic  ray  theory,  ll^e  reference  point  in  the  evaluation  of 
these  changes  is  a  medium  where  there  is  no  temperature  gradient.  For  vertical 
rays  where  Officer's  eqn  (2)  is  invalid,  analysis  of  e<p  (4)  iidicates  no  intensity  or 

phase  change.  Equation  (4)  also  fails  at  6  =  ^  which  corresponds  to  a  turning 

point.  Ray  solutions  always  fail  at  turning  points.  Since  the  acoustic  velocity  must 
be  positive,  for  0g  not  equal  to  0,  the  horizontal  range  L  cannot  be  equal  to  0.  The 
expressions  obtaiired,  eqns  (1, 3, 5, 6),  enable  one  to  obtain  intensity  and  phase 
changes  due  to  variations  in  the  mean  temperature  gradient,  in  regions  where  the 
ray  approximation  is  valid. 

Section  4  and  beyond  discusses  models  to  obtain  acoustic  estimates  due  to  the 
random  part  of  the  temperature.  Equation  (8)  is  the  expression  for  the  scattered 
pressure,  using  a  stochastic  wave  model.  Numerical  calculatioiu  of  the  angular 
distribution  of  the  scattered  pressure  for  high  frequency  sound  (Figs.7, 8)  both  for 
anisotropic  and  isotropic  media  indicate  scattering  being  confined  to  an  angle 
close  to  the  direction  of  propagation.  For  attenuation  in  an  isotropic  medium,  an 
expression  has  been  available  for  some  time,  eqn  (12)  [Chernov,  1960,  p55].  In  an 
anisotropic  medium,  attenuation  depends  on  the  direction  of  propagation.  When 
the  direction  of  propagation  is  along  the  Z  axis,  attenuation  is  given  by  eqn  (13). 
At  1  MHz  the  attenuation  is  of  the  order  of  5X10"  ^  dB/ m.  The  magnitudes  for 
the  isotropic  medium  are  of  the  same  order.  Attenuation  varied  as  kl  for 
frequencies  above  100  kHz. 

Section  4.2  examines  analytically  amplitude  and  phase  fluctuations  at  a  receiver 
separated  by  a  horizontal  range  L  from  a  continuous  monochromatic  acoustic 
source.  The  study  is  in  the  unsaturated  region  of  the  A  -  4>  space  and  uses 
stochastic  wave  and  ray  models.  For  an  isotropic  medium,  with  very  much 
greater  than  1  and  range  L  very  much  greater  than  a^,  expressions  (18)  and  (19) 
for  acoustic  fluctuations  have  been  available  in  the  wave  model  [Chernov,  1960, 
p83].  Equations  (14-17)  are  the  new  results  for  acoustic  fluctuations  in  an  isotropic 
medium,  valid  at  all  ranges  and  frequencies.  Numerical  comparisons  between  the 
two  models  show  (see  Fig.  9)  that  the  Chernov  model  is  valid  only  for  L  very 
much  greater  than  a^,  and  k^a^  very  much  greater  than  1.  Equations  (14, 15, 27- 
30)  give  expressions  for  acoustic  fluctuations  in  an  anisotropic  medium,  which  are 
valid  at  all  ranges,  and  frequencies.  The  anisotropic  case  has  not  been  treated 
previously. 

Expressions  (24)  and  (25)  for  acoustic  fluctuations  in  the  stochastic  ray  model 
have  been  available  only  for  an  isotropic  medium  and  at  large  range,  where  L  is 
very  much  greater  than  [Chernov,  I960,  p29, 34].  Equations  (22, 23)  are  the 
new  results  for  the  acoustic  fluctuations  in  the  ray  approximation  valid  at  all 
ranges  for  an  isotropic  medium.  Numericai  comparisons  between  the  two  ray 
models  (Fig.  10)  show  the  Chernov  mode)  is  valid  only  for  L  very  much  greater 
than  0^ .  Equation  (22)  with  replacing  gives  the  expression  for  the  rms 
phase  fluctuation  in  the  stochastic  ray  model  for  an  anisotropic  medium.  The 
anisotropic  case  has  not  been  treated  previously. 
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Sections  4.3  amd  4.4  presents  numf^rical  estimates  of  intensity  and  phase 
fluctuations.  Calculations  at  100  kHz,  wich  L  equal  to  100  m  and  a  temperature 
fluctuation  of  7  X 10“^  C  deg,  show  the  rms  relative  intensitv  fluctuations  (dB)  and 
rms  phase(rads)  for  an  anisotropic  medium  are  on  the  average  Irigber  by  a  factor 
of  5-10,  and  3-6  respectively,  when  compared  with  an  isotropic  medium,  where 
the  vertical  correlation  lengths  are  the  same,  and  the  radial  correlation  length  is 
about  50  times  larger.  For  these  parameters  the  rms  phase  and  the  rms  relative 
intensity  for  an  isotropic  medium  are  of  the  order  of  0.1  radians  and  0.3  dB 

respectively.  1  he  variation  of  the  rms  relative  intensity  and  rms  phase  with  , 
which  is  dependent  on  the  magnitude  of  the  temperature  fluctuation,  is  linear 
both  for  the  isotropic  and  anisotropic  media.  The  variation  of  the  acoustic 
fluctuations  as  a  function  of  frequency  for  large  obtained  using  the  wave 
model,  show  that  the  agreement  of  the  rms  phase  with  the  ray  result  occurs  at  a 
lower  frequency  than  for  the  relative  intensity.  A  shorter  range  gives  better 
agreement  with  the  ray  model,  since  an  increase  in  range  increases  the  diffraction 

parameter  and  .  Comparison  with  the  statistical  ray  model  shows,  that  in  the 
isotropic  case  diffraction  effects  become  important  at  correlation  lengths  that  are 

less  than  approximately  • 


Finally,  section  5  examines  the  evaluation  of  theoretical  estimates  for  phase 
structure  function,  using  the  ray  model.  Agreement  in  the  phase  fluctuation 
estimates  between  the  wave  and  ray  models  is  within  5%  for  an  isotropic  medium 
with  L  equal  to  100  m,  if  k^Pf^  is  greater  than  500.  The  same  agreement  occurs  in 
the  anisotropic  medium  if  and  A:fjp„  is  greater  than  600.  At  a  range  of  3  m 
and  a  frequency  of  1  MHz,  the  results  of  the  phase  fluctuations  from  the  two 
models  show  that  they  are  indistinguishable  graphically  (Fig.  18).  Equations  (39) 
cind  (40)  are  the  expressions  of  the  phase  structure  function  for  an  isotropic  and 
anisotropic  medium  respectively,  valid  at  ranges  L  »  and  at  separations 
2  =  In  other  words  eqns  (39),  (40)  are  expressions  for  the  phase  structure 
function  at  large  range  and  for  small  separations  between  receivers.  One  can  also 
obtain  eqn  (39)  in  the  ray  approximation,  using  the  spatial  auto  correlation 
function  of  the  phase  fluctuations  [Chernov,  1960,  plO?].  Equations  (34,  35)  and 
eqns  (36,  37)  are  the  expressions  for  the  PSF  of  isotropic  and  anisotropic  media 
respectively,  valid  at  short  range,  and  at  all  separations  between  receivers.  At 


short  range,  the  limiting  value  of  ^  for  large 


«t>f 


is  1  for  the  anisotropic 


medium,  and  2  for  an  isotropic  medium.  At  large  range  where  L  is  very  much 
greater  than  P(,,  both  approach  the  limiting  value  of  2. 
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Appendix 

List  of  Symbols 


f 

/ 

fs 

c 

k 

t 

T 

S 

r 

z 

S 

sm0o 
P  =  -^ 

s 

dJ 

9 

<& 

y 

M 

K 

A 

4' 


frequency 
inertial  frequency 
Brunt  Vaisald  frequency 
sound  velocity 

wave  number  = 

c 

tine 

degree  Centigrade 
salinity  in  ppt 
positional  coordinate 
depth  in  m 
velocity  gradient 

Cq  speed  of  sound  near  source;  By  ray  angle  at  source 

ray  path 
Intensity  loss 
phase 

root  mean  square  phase  fluctuation  in  the  ray 
approximation 

used  in  equation  (6) 

mean  square  scattered  pressure 
diffe'^nt  temperature  profiles 
wavelength  of  sound 
diffraction  parameter 

mean  square  fractional  sound  speed  fluctauation 
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L 


N,{r-r') 

Nir  -r'  ) 


amplitude  of  received,  transmitted  sound  wave 


covariance 


correlation  coefficient 


K  —  k—  ;  /'  V  denote  direction  of  incident  and  scattered  sound 


F(K) 


Fourier  transform  oi  N(r-r') 

frequency  and  angular  part  of  the  mean  scattered 
pressure 

total  attenuation  in  dB/m 
vertical  and  radial  correlation  length 
isotropic  correlation  length 

mean  square  relative  log  amplitude  fluctuation 

mean  square  phase  fluctuation 

integrals  that  appear  in  fluctuation  formulae 


IA,IB,IC,ID  defined  by  equations  (29)  and  (30) 

S(g),C(tj)  Fresnel  sine  and  cosine  integrals 

Eq  energy  emitted  per  unit  time  per  unit  solid  angle 

PSF(x)  phase  structure  function 

1,^22  >^33  defined  by  equations  (35),  (37)  (38) 
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